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ABSTRACT 

New observations over the next few years of the emission of distant objects will help 
unfold the chapter in cosmic history around the era of the first galaxies. These obser- 
vations will use the neutral hydrogen emission or absorption at a wavelength of 21-cm 
as a detector of the hydrogen abundance. We predict the signature on the 21-cm signal 
of the early generations of galaxies. We calculate the 21-cm power spectrum including 
two physical effects that were neglected in previous calculations. The first is the redis- 
tribution of the UV photons from the first galaxies due to their scattering off of the 
neutral hydrogen, which results in an enhancement of the 21-cm signal. The second is 
the presence of an ionized hydrogen bubble near each source, which produces a cutoff 
at observable scales. We show that the resulting clear signature in the 21-cm power 
spectrum can be used to detect and study the population of galaxies that formed just 
200 million years after the Big Bang. 

Key words: galaxies:high-redshift - cosmology:theory 



1 INTRODUCTION 

An important milestone in the evolution of the Universe is 
the appearance of the first luminous objects, which begins 
the era of heating and ionization of the intergalactic medium 
(IGM) , referred to as cosmic reionization. A promising probe 
of the cosmic gas up to the end of reionization is the hyper- 
fine spin-flip transition of neutral hydrogen (H I) at a wave- 
length of 21-cm. Observations of the redshifted 21-cm power 
spectrum as a function of wavelength and angular position 
can provide a three-dimensional map of t he H I distribution 
l|Hogan fe Reeslll979l ; iMadau et al.lll997l) . 

The 21-cm line of the hydrogen atom results from the 
transition between the triplet and singlet hyperfine lev- 
els of the ground state. Calculations of 21-cm absorption 
begin with the spin temperature T s , defined through the 
ratio between the number densities of hydrogen atoms, 
Tii /no = 3 exp(— T+/T s ), where subscripts 1 and corre- 
spond to the excited and ground state levels of the 21-cm 
transition and T* = 0.0682 K corresponds to the energy dif- 
ference between the levels. The 21-cm spin temperature is on 
the one hand radiatively coupled to the cosmic microwave 
background (CMB) temperature T T , and on the other hand 
coupled to the kinetic gas t emperature Tk through collisions 
(All ison fc Dalgarnol Il969f ) or the absorption of Lya pho- 
tons, as discussed b elow. For the concor dance set of cosmo- 
logical parameters (S pergel et alj 120071 ) , the mean bright- 
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ness temperature on the sky at redshift z (relative to the 
CMB itself) is T b = 28 mK y/(l + z)/W [(T s - T 7 )/T s ] xm 
l|Madau et al.lfl997l ) where xhi is the mean neutral fraction 
of hydrogen. 

After cosmic recombination, down to z ~ 200 the re- 
maining free electrons kept Tk close to T 7 . Afterward, the gas 
cooled adiabatically and atomic collisions kept the spin tem- 
perature T s coupled to Tk, both of them lower than T 7 , thus 
creating a 21-cm signal in absorption. As the universe con- 
tinued to expand and rarefy the gas, the radiative coupling 
between T 7 and T s started to dom inate again, and the 21- 
cm si gnal faded. Starting at z ~ 66 (|Naoz. Noter fc Barkanal 
2006), the UV photons from the first galaxies emitted be- 
tween the Lya and Lyman limit wavelengths propagated 
freely in the IGM. They redshifted or scattered into the Lya 
resonance, thus coupling T„ to Tk through the WF effect 
l|Wouthuvsen|[l952l ; [Field 1958) whereby atoms re-emitting 
Lya photons can de-excite into either of the hyperfine states. 

Fluctuations in the Lya radiation emitted by the 
first galaxies led to fluctuations in the 21-cm signal 
IjBarkana fc Loebl2005bl ). until the WF effect saturated (i.e., 
finished dropping T s to Tk). In particular, the biased fluc- 
tuations in the number density of galaxies, combined with 
Poisson fluctuations in the number of galaxies and the 
1/r 2 dependence of the flux, caused large fluctuations in 
the Lya background. This Lya background is composed of 
two parts, the photons that directly redshift into the Lya 
resonance, and those produced during the atomic cascade 
from higher Lyman series photons (Ly7 and above). About 
30% of the total number of photons emitted from the first 
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stars between Lyy and the Lyman limit are co nverted to 
Lya i|Pritchard fc Furlanettol 120061 ; iHiratal I2006T ). The 21- 
cm power spectrum that arises from the fluctuations in 
the Lya background can be used to pr obe the first sources 
of light and their effect on th e IGM dMadau et al.lll997 ; 
Chen fc Miralda-Escudel 12004 IBarkana fc Loebl l2005al lb; 



Chuzhov. Alvarez fc Shapiroll2006i r 



In this Letter we focus on two effects that have not 
previously been included in calculations of statistics of the 
21-cm brightness temperature; these effects dramatically in- 
fluence the Lya flux fluctuations and therefore the 21-cm 
power spectrum. The first is the diffusion of photons through 
the neutral hydrogen gas. Since the optical depth of Lya is 
around a million, photons scattered on the surrounding gas 
cause each so urce to appear as a halo in the sky rather than a 
point source i|Loeb fc Rvbickilll999T ). In the absence of scat- 
tering a photon will redshift to Lya at some distance from 
the source, but in the presence of scattering the emitted pho- 
tons scatter ba ck and forth and thus re a ch Lya much closer 
to the source JChuzhov fc Zhendl2007l : ISemelin fc Combel 
120071 ). Thus, the Lya flux that a random point in the uni- 
verse observes is dominated by nearby sources. As a result, 
the fluctuations in the Lya flux are larger and enhance the 
21-cm power spectrum. The second effect is caused by the 
presence of ionized hyd rogen (H II) around each source (e.g., 
IBarkana fc Loebl l200ll ). In this region photons can redshift 
below the Lya resonance without encountering H I, and thus 
not participate in the WF effect. Thus, a typical H atom 
(which by definition is not within an H II region) receives 
no Lya photons from sources closer than the typical size 
of an H II bubble. Since galaxy number fluctuations below 
this scale do not contribute to Lya fluctuations, we expect 
a cutoff in the 21-cm power spectrum around this scale. 

Our calculati ons are made in a A CDM universe match- 
ing observations (|Spergel et al.ll2007l ). with a power spec- 
trum normalization erg = 0.826, Hubble constant Ho = 
68.7 km s _1 Mpc -1 , and density parameters Q, m — 0.299, 
SIa = 0.701, and Slf, = 0.0478 for matter, cosmological con- 
stant, and baryons, respectively. 



2 LYMAN SERIES SCATTERING 

We assume a uniform neutral IGM that ex hibits pure Hub- 
ble ex pansion around a steady point source (|Loeb fc Rvbickil 
1999). The surrounding IGM, mainly H I, is optically thick 
to the UV radiation that the sou rce emits. We use a Mo nte- 
Carlo method for the scattering l|Loeb fc Rv bicki 199§), for 
photons of various initial frequencies between Lya and the 
Lyman limit. The scattering redistributes the distances at 
which the photons transform to Lya, affecting strongly the 
photons originally emitted between Lya and Ly/3 . Due to 
flux conservation, the enhancement of the flux at small scales 
(by up to a factor of ~ 5 compared to the optically thin case) 
is balanced by a steep drop on large scales (where v — >Ly/3). 
Lya photons produced via cascading from photons emitted 
above Ly7 are much less affected (see below), but we fully 
include these cascades assuming tha t the resulting Lya pho- 
tons are emitted at the line center (|Furlanetto fc Pritchardl 
2006). We account for several features that substantially af- 
fect the redistribution of photons compared to previous cal- 
culations that gave a divergent r _1,/3 enhancement near the 



source l|Chuzhov fc Zhendl2007l ; | Semelin fc Combes! I2007T I. 
First, we explicitly include the H II region around each 
source, which results in enhanced Lya scattering just out- 
side the H II region; the enhancement has a shallower rise 
to a peak value and then drops to zero right near the H II 
region, because of the loss of photons that are scattered back 
into the H II region. Second, each hydrogen atom receives 
Lya photons from sources with an effective spectrum that 
drops sharply at frequencies approaching Ly/3; photons at 
such frequencies are emitted by time-retarded sources that 
formed in a universe with fewer galaxies, effectively reduc- 
ing the emission rate (see also section 4 below). We find 
that scattering strongly enhances this drop in the spectrum. 
Finally, we also account for the variation of the Hubble con- 
stant with redshift, which substantially affects the photons 

emitted at v — >Ly/3. 

Following lLoeb fc Rvbickil (ll999T ) we define a dimen- 
sionless variable f = r/r*, where r* is the radius for which 
redshifting a Lya photon due to Hubble expansion would 
produce a frequency with a remaining optical depth of unity 
(from there out to infinity, in the direction away from the 
line center). For example, using the standard cosmological 
parameters, r* ~ 1 Mpc (physi cal), independent of red shift 
(for a source at high redshift: lLoeb fc Rvbickil <|l999l )). In 
Figure [T] we show the ratio between the flux of the source 
as a result of our calculation as described above and the 
flux for the case with no redistribution due to scattering, 
i.e., where the photons simply redshift until they are ab- 
sorbed at line center. For photons emitted between Lya and 
Ly/3 we have plotted the ratio as a function of the vari- 
able x = (f - fmi)/(r m ax ~ rim), where r max is the radius 
at which an emitted Ly/3 photon redshifts all the way to 
Lya, and the various radii are expressed in units of r*. At 
z = 20, r* = 22.8 comoving Mpc and r max = 321.4 co- 
moving Mpc. Showing the results as a function of the vari- 
able x provides a better view of the behavior of the flux 
very near to the H II region^. We have considered in this 
Figure a Population II spectrum source, i.e., we approxi- 
mated the emissivity within this fr equency region as ~ v as , 
where for Population II a s = 0.14 |Barkana fc Loeb| [2005b: 
iBromm. Kudritzki. fc Lo cb 2001). We have also performed 
the calculation for a spectrum appropriate for Population 
III stars (see below). For photons cascading from higher- 
order Lyman lines, the atomic constants would make the 
enhancement occur much closer to the H II region edge; 
e.g., even for Ly7, r* = 0.152 comoving Mpc, r max = 40.3 
comoving Mpc, and a substantial enhancement occurs only 
at x 10 -5 (with the corresponding definition of x for Ly7), 
well within the thickness of the partial ionization shell of the 
H II region, where the flux is in practice strongly suppressed. 



3 21CM FLUCTUATIONS 

In general, fluctuations in Tf, can be sourced by fluctuations 
in the gas density, temperature, neutral fraction, radial ve- 



1 Note that the numerical noise seen in the Figure very near the 
H II region occurs where the flux is not significantly enhanced, 
and thus does not affect our results. Furthermore, this region is 
within the thickness of the shell of the H II region, and thus is 
even less important in practice. 
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Figure 1. The ratio between the flux as a result of our calculation 
as described in the text and the flux for the case with no scatter- 
ing redistribution, for photons emitted between Lya and Ly/3 at 
z = 20. We consider an H II region of size 0.023, 0.23, and 2.3 co- 
moving Mpc (open squares, triangles, and crosses, respectively). 
Assuming Population II stars, the filled squares (light green) are 
for the case of a simple power-law spectrum (with rjjn = 0.023 
Mpc), while the other curves assume the correct effective spec- 
trum with time-retarded sources. For the effective spectrum, the 
flux is low at large distances, so the scattering (which conserves 
the total photon number) must lower this flux especially strongly 
in order to compensate for the enhancement at small distances. 



locity gradient, and Lya flux from galaxies. To linear order, 
perturbing the equation of the brightness temperature, we 
have 8Tb = PbS b + PtSt + 5 X + f3 a S a — 5d rVr , where each 
Si describes the fractional perturbation on the quantity i; b 
denotes the baryon density, T the gas temperature, x the 
neutral fraction, a the Lya coefficient x a (which is simply 
the Lya flux measured in units so that Lya flux fluctua- 
tions give the largest 21-cm fluctuations when x a — 1), 
and d r v r the line of sight velocity gradient [which creates 
an apparent density fl uctuation l|Bharadwai fc AlH |2000| ; 
iBarkana fc Loebl l2005al )]. The associated time-dependent 
Pi coefficients are /3 6 = 1 + [x c /x to t}, Pa = x a /x to t, 
Pt = [T 7 /(T fe -T 7 )] + (x c /xtot)*[rflog«;i-o/rflogT fc ], where 
itot = 2tot(l + ztot), xtot = x a + x c , and x c and ki_q de- 
scrib e the collisional 21-cm excitation ( Allison fc Dalgarnol 
1 19691 ; IBarkana fc Loebl l2005al ; iFurlanetto et alj I2006T ). We 
consider the high-redshift regime before significant cosmic 
reionization or stellar heating. Due to Compton scattering 
of the CMB photons with the remaining free electrons after 
cosmic recombination, the baryon density and temperatu re 
fluctuations are not proportional (|Naoz fc B arkana 2005). 

A Fourier transform yields the brightness temperature 
fluctuation in k- space (where S denotes the transform of 
each 3 quantity) l|Barkana fc Loebl l2005al ): 



8T b {k,t) = (fi 2 r lM + P) S b (k) + pj a {k) 



(1) 



where the ratio 



5:i](^' 



[(d/dt)S b ]/(HS b ), P(k,z) 



Pb + r[ T: s]Pr in terms of ry T .^ = Sr/Sb, and n = cos8k 
in terms of the angle 6t of the wavevector k with respect to 
the line of sight. 

We denote by Ps b (k) and P a {k) the power spectra of the 
fluctuations in baryon density and in the Lya flux, respec- 
tively, and the power spectrum Ps b - a as t he Fourier trans- 
form of their cross-correlation function l|Barkana fc Loebl 
2005b). The 21-cm power spectrum can then be written in 
the form P Tb (k) = fi 4 P^(k) + M 2 i>(*0 + *>(fc) , where 
each of the three [i coefficients ca n be separately measure d 
from their different jx dependence ( Barkana fc Loebll2005al ). 
This angular separation of power makes it possible to detect 
separately different physical aspects influencing the 21-cm 
signal. In terms of our definitions above, the coefficients are 



P A k ) = r [h . A P Sb {k) 

PA k ) = 2r lS:5] {pP Sb (k) + PaP6 b ^a(k)] 
P M o(fc) = p 2 P Sb {k) +p 2 a P a {k) +2PPaP Sb -a(k) 



(2) 



It is possible to probe whether some sources of P a are un- 
corrected with Sb\ sources that are linear functionals of the 
baryon densit y distribution do not co ntribute to the follow- 
ing quantity ( Barkana fc Loe b 2005a): 



,(*) 



4P„ 



= Pl 



p2 



(3) 



4 THE Lya FLUX OF GALAXIES 

We summarize here t he calculation of the Lya flux 
IIBarkana fc Loebll2005bl ). The intensity of Lya photons ob- 
served at a given redshift is due to photons that were orig- 
inally emitted between the rest-frame wavelengths of Lya 
and the Lyman limit. Photons that were emitted below Ly/3 
by a source simply redshift until they are seen by an atom at 
z at the Lya wavelength. Such photons can only be seen out 
to a distance corresponding to the redshift z max (2), where 
[1 + 2 m ax(2)]/(l + z) = (32/27), the ratio of Lya to Ly/3 
wavelengths. Photons above the Ly/3 energy redshift until 
they reach the nearest atomic level n. The neutral IGM is 
opaque even to the higher levels and so the photon is repeat- 
edly scattered, with a ~ 30% chance of being downgraded 
to a Lya photon and continuing to scatter, except that the 
chance is zero fo r a Ly/3 photon (jPritchard fc Furlanettol 
l2006l ; lHiratal200t3 ). To be seen at the Lya resonance at z, the 
photon must have been emitte d below a maximal red shift 
2 m „(?i). The intensity is then (Barkana & Locb 2005b): 



J a — 



(1 + ^) 2 



47T 

n max 



(4) 



frecyc (^) 



cdz 



H(z 



jre(y n ,z )/ scat (n,r) 



where absorption at level n at redshift z corresponds to an 
emitted frequency v' n at z', and e i s the photon emissivity. 
We have included the factor / reC yc (jPritchard fc Furlanettol 
l2006l : lHiratall2006l ) which is the fraction of photons absorbed 
at level n that are eventually recycled into Lya photons in 
their subsequent cascade. The new correction factor / scat 
is the overall factor by which the Lya flux is changed due 
to the scattering-induced redistribution of photons that are 
emitted with frequencies between levels n and n + 1 from 
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sources at a comoving distance r from the final destination 
at redshift z (where r is a function of z and z'). Also new to 
calculations of the 21-cm power spectrum is the lower limit 
of the integral (zhii), which accounts for the H II regions; 
in order to receive Lya photons from sources closer than the 
typical size of an H II region, an H atom would have to be 
inside an H II bubble - a contradiction. 



5 THE 21-CM POWER SPECTRUM 

As mentioned above, the re are two separate sour ces of fluc- 
tuations in the Lya flux ([Barkana fc Lo cb 2005b]). The first 
is density inhomogeneities. Since gravitational instability 
proceeds faster in overdense regions, the biased distribution 
of rare galactic halos fluctuates much more than the globa l 
dark matter density ((Kaiser 1987; Barkan a fc Loebll2004h . 
When the number of sources seen by each atom is relatively 
small, Poisson fluctuations provide a second source of fluctu- 
ations, statistically independent of the first source to linear 
order. Unlike typical Poisson noise, these fluctuations are 
correlated between gas elements at different places, since 
two nearby elements s ee many of the same sou rces, though 
at different distances (jBarkana fc Loebl l2005bT l. Note that 
although each hydrogen atom receives some Lya flux from 
sources as far away as 300 comoving Mpc, the fluctuations 
in flux are actually relatively large because a significant por- 
tion of the flux comes from nearby sources. 

Due to the geometrical dependence of the flux, the per- 
turbation in the Lya flux due to biased density fluctuations 
is a linear functional of the underlying density fluctuation 
field, i.e., the resulting contribution to <5 Q (k) is related to the 
Fourier transform of the total density perturbation 5 to t(k) 
by multiplication by an effective window function W(k) 
(Bark ana fc Loe b 2005b). The total density perturbation is 
defined as the mass-weighted mean of the dark matter and 
baryon perturbations; we also define r[S:S b ] = Stot/Sb- Under 
the conditions considered in this paper, the three observable 
power spectra (eq. [2} can be used to study separately the 
two sources of Lya fluctuations. In particular, the quantity 

P flux _,(fc) = P^ 2 -^Lp^ = 2r^. s] (3 a r [S:Sb] W(k) P Sb (k)(5) 

r [5:S] 

is proportional to the biased-density contribution to P a , 
while P un -s(k) equals the Poisson contribution to P a ex- 
cept for a factor of /3„. 



6 OBSERVABLE SIGNATURE 

The coupling of the spin temperature to the kinetic temper- 
ature through the WF effect requires a relatively low cosmic 
Lya background, and is exp ected to occur well before the 
end of cosmic reionization dMadau et alj ll997T l . We illus- 
trate our results assuming that the Lya coupling transition 
achieves x a = 1 at z — 20; this requirement determines 
the star formation efficiency. We assume that galaxies form 
within all dark matter halos above some minimum mass (or, 
equivalently, a minimum circular velocity V c = yj GM/R in 
terms of the virial radius R). We consider a wide range of 
possible values: galaxies that form through molecular hy- 
drogen cooling (V c = 4.5 km/s), atomic cooling (V c = 16.5 



km/s), or a minimum mass ten times larger (V c = 35.5 
km/s) due to feedback effects in low-mass halos. We con- 
sider a single H II region size (calculated as having the flux- 
weighted mean volume); this size is affected by the fraction 
/ csc of stellar ionizing photons that escape into the IGM, 
and by the stellar population. We consider two extremes, 
a stellar initial mass function as observed locally (Pop II) 
or that expected for the very first stars (Pop III; 100 solar 
mass, zero-metallicity stars). 

In Figure [2] we show in two representative cases the 
two separately-measurable power spectrum terms that are 
due to Lya flux fluctuations arising from biased density 
fluctuations (-Pfl ux -a) or from Poisson fluctuations (P un -s)- 
Each term (Pfl ux _,5 or P un _,5) starts small on large scales 
and rises with k, forming a peak before dropping (and 
then oscillating) on the scale of the H II region. Scatter- 
ing modifies P U n-s more substantially than Paux-s, since 
the former is more strongly dominated by fluctuations on 
small scales. Note that our assumption of a single H II re- 
gion size is reasonable since galaxies at these redshifts are 
rare and almost all lie in halos close to the minimum halo 
mass. Also, sufficiently early on in reionization, galaxy clus- 
tering plays a relatively minor role in determining bubble 
sizes. In particular, in the cases shown in Figure [2] the 
cosmic ionized fraction is ~ 10~ 3 , and bubbles of sizes be- 
tween the minimum size (that of a single-source bubble) and 
twice that size contain most of the ionized volume [2/3 for 
V c = 16.5 km/s, and 1/2 for V c = 35.5 k m/s, using the 
iFurlanetto. Zaldarriaga fc Hernquistl {2004) model]. While 
the small-scale ringing seen in Figure [2] may be smoothed 
out by a mild scatter in H II region sizes, the overall shape 
and in particular the peak of each power spectrum term are 
robustly determined. 

We summarize the main features of the two power 
spectra through the peak positions (fcp"*^ 15 and fc™^ 5 ) 
and heights (T^~ s and T^ s ). Now, since the flux term 
would have a peak eve n without the H II region cutoff 
l|Barkana fc Loe b 2005b), the characteristics of its peak are 
fairly insensitive to the H II region size; thus, we also con- 
sider the position of its cutoff, specifically the lowest k value 
(above k^^ 1 ) where the power spectrum drops to 1 mK: 

fc imK S - Both fc pca~k 4 and essentially measure the size 

J?hii of the H II region (figure [3]). In particular, the product 
fe p^k * R Bu ~ 0.6 and kf^~ s * R ml « 1-6, each to within 
a factor of 1.5 over a range of three orders of magnitude 
of -Rhii- On the other hand, the position and height of the 
peak of Pflux-6 are relatively insensitive to -Rhii and thus 
observing them would constitute a clear consistency check 
with the theory. The Poisson peak height T^~^ s measures 
the average number density of galaxies and thus depends 
mainly on the minimum mass of galactic halos. 

We note that if significant X-ray heating happened to 
occur simultaneously with the Lya coupling transition, then 
similar fluctuations in the X-ray flux from galaxies would 
generate 21-cm fluctuations th at are even somewhat larger 
(jPritchard fc Furlanetto~ll2007h : while not included in pre- 
vious studies, the X-ray fluctuations should similarly be en- 
hanced by scattering and should also show a cutoff at the 
H II bubble scale. While X-ray heating could possibly oc- 
cur early, Lya heating i s insignificant except after the Lya 
coupling fully saturates (|Chen fc Miralda-Escudel [20041 ). 
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Figure 2. 21-cm power spectrum P (shown in terms of the 
brightness temperature fluctuation) as a function of the comoving 
wavenumber k. We consider Pfl ux _f (top panels) and P un _$ (bot- 
tom panels). We compare the previous result ( Barkana fc Loebl 
l2005al ; iPritchard fc Furlanettoll2006l; lHiratall200rJ) (short-dashed 
curves) and the result corr ected to use the precis e density and 
temperature power spectra llNaoz &: B arkana 2005) (long-dashed 
curves) to the full calculation with the H II region cutoff (solid 
curves, with negative portions dotted and shown in absolute 
value; the higher curve of each pair also includes the redistri- 
bution of photons due to scattering). We consider galactic halos 
with a minimum circular velocity V c = 16.5 km/s (left panels) or 
35.5 km/s (right panels), assuming Pop II stars with / csc = 0.3. 

The five quantitative characteristics we have focused 
on typically occur at relatively large scales - a fraction of 
a Mpc, or ~ 10 arcseconds at z = 20 - and correspond to 
relatively large (but still linear) fluctuations: 1-10 mK on a 
mean background of Tt ~ —100 mK. Thus, the predictions 
are theoretically robust, and require observational capabili- 
ties only somewhat beyond those of the radio arrays being 
currently constructed. Given current capabilities, 21-cm cos- 
mology is the most promising method for firmly detecting 
and studying the properties of some of the earliest galaxies 
that ever formed. 
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Figure 3. Main features of 21-cm power spectra from Lyct flux 
fluctuations as a function of the comoving H II region size. We 
show (top panel) k^^~ S (dashed curves), k^~^ (solid curves), 

and fefj^K * (dotted curves), as well as (bottom panel) T^^T S 
(dashed curves) and (dotted curves); see text for defini- 

tions. These curves are the result of a sweep of the parameter 
space, made in six segments each of which varies / esc from 1— 
100% thus covering a range of H II region sizes shown by one of 
the horizontal bars (bottom panel): V c = 4.5 km/s with Pop II 
stars or Pop III stars, V c = 16.5 km/s with Pop II or III, and 
V c = 35.5 km/s with Pop II or III stars, from left to right. 
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